Point symmetries in the Hartree-Fock approach: Symmetry-breaking schemes 
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We analyze breaking of symmetries that belong to the double point group D™ (three mutually 
perpendicular symmetry axes of the second order, inversion, and time reversal). Subgroup structure 
of the D™ group indicates that there can be as much as 28 physically different, broken-symmetry 



mean-field schemes - 



starting with solutions obeying all the symmetries of the D™ g rou P> through 



break all of the D™ symmetries. 



26 generic schemes in which only a non-trivial subgroup of D 2h is conserved, down to solutions that 

Choices of single-particle bases and the corresponding structures 
of single-particle hermitian operators are discussed for several subgroups of D™ . 

PACS numbers: 21.60.-n, 21.60.Jz, 21.10.Ky 



I. INTRODUCTION 

One of the salient features of the mean-field approach 
to many-fermion (e.g., nuclear) systems is the sponta- 
neous symmetry breaking. The symmetry of a mean- field 
state is called broken, if the solution of the Hartree-Fock 
(HF) or Hartree-Fock-Bogolyubov (HFB) self-consistent 
equations do not obey symmetries of the original many- 
body Hamiltonian |IJ . This happens when the calculated 
mean-field energy of the system is lower for states which 
break a symmetry than that for unbroken symmetries. 
Such a mechanism depends on the physical situation and 
is governed by the Jahn- Teller effect Without going 
into details, let us recall that the spontaneous breaking of 
an original symmetry is usually accompanied by a signif- 
icant decrease in the single-particle level density at the 
Fermi energy. Hence, the doubly magic nuclei can be 
safely described by imposing conservation of the spheri- 
cal symmetry, while this symmetry should be allowed to 
be broken in the open-shell systems. 

One of the simplest examples in this context is that of 
the breaking of the translational symmetry. The related 
mechanism is present, e.g., in the nuclear shell-model. 
Indeed, within the framework of the shell model, inter- 
acting nucleons are assumed to move in a common mean 
field that is localized in space and consequently they can- 
not be described by eigenstates of the momentum opera- 
tor (plane waves). In other words, the wave functions of a 
nucleus cannot be approximated by uncorrelated single- 
particle plane waves - this can only be attempted for an 
infinite system, i.e., for the nuclear matter. The use of a 
shell-model, space-localized wave function simply reflects 
the correlations present in the system. In this example, 
the correlations ensure that it is improbable to find two 
nucleons of a nucleus at large relative distances apart. 

In nuclear structure physics one can easily identify the 
use of various broken symmetries in a description of well- 
defined, observable effects. For instance the rotational, 
parity, time-reversal, and gauge symmetry breaking were 
introduced to describe the deformations, octupole corre- 
lations, nuclear rotation and pair correlations, and com- 



bined effects thereof. At present, we approach the situa- 
tion where the mean-field calculations can be performed 
without explicitly using any of the mean-field symme- 
tries. Several such approaches have already been imple- 
mented 0-0], although very few calculations for specific 
physical problems have been done to date. 

One could, in principle, perform the mean-field calcu- 
lations without assuming a priori any symmetry, and let 
the dynamics choose those discrete symmetries which are, 
in a specific situation, broken, and those which remain 
obeyed. Obviously, by choosing such an approach we can- 
not profit from simplifications possible when it is known 
beforehand that some symmetries are obeyed/disobeyed. 
However, following the general guidelines provided by 
the Jahn- Teller mechanism one usually can make a rea- 
sonable choice of obeyed/disobeyed symmetries. Such a 
choice is dictated by the properties of the many-body 
Hamiltonian and by the classes of phenomena which one 
wants to describe - it usually facilitates the calculations 
markedly. In all those cases the analysis presented in this 
article provides us with the mathematical means for con- 
structing the algorithms optimally adapted to the sym- 
metries of the problem in question. 

In the preceding article ||, we have presented prop- 
erties of the single point group and double point 
group D™, that can be built from operators related to 
the three mutually perpendicular symmetry axes of the 
second order, inversion, and time reversal. We have also 
discussed their roles in the description of even and odd 
fermion systems, respectively, their representations, and 
the symmetry conditions induced by the conserved Dj h 
or D™ symmetries on the local densities and electro- 
magnetic moments. 

By considering the D™ double point group we focus 
on quantum objects that are in general non-spherical, but 
can have one or more symmetry axes and/or symmetry 
planes. Obviously, any nuclear many-body Hamiltonian 
of an isolated system is time-even and rotationally in- 
variant. In the present paper we do not aim at analyzing 
the conditions under which these symmetries are broken 
spontaneously, with one or another symmetry element of 
the D™ group still being conserved in the HF solution. 
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Instead, we present a classification of all such possibili- 
ties, and discuss the resulting properties of the mean-field 
Hamiltonians and single-particle wave functions. For a 
review of applications of point symmetries to a descrip- 
tion of rotating nuclei see the recent study in Ref. 0. 

Our goal is thus twofold: First, in Sec. |l| we discuss 
all possible physically meaningful subgroups of D™ , and 
classify the corresponding physical situations from the 
view-point of the conserved D™ symmetries. In many 
applications to date, specific choices of conserved and 
broken D™ symmetries have been made [p[-^4[, how- 
ever, here we aim at a complete description of all achiev- 



Dj h : {E,P,R fc ,S fc ,f,P T ,R^,Sn, 



(3) 



able symmetry-breaking schemes. Second, in Sec. [II 
we review and discuss practical aspects of structure of 
the mean-field operators under specific D™ group op- 
erations. This essential question has been explicitly or 
tacitly addressed in most approaches using the deformed 
mean-field theory; our aim here is to present exhaustive 
list of options pertaining to all the D™ symmetry con- 



ditions. Finally, conclusions are presented in Sec. IV 



II. SUBGROUPS OF Dj h AND D^P AND THE 
SYMMETRY BREAKING 

The single group Dj h and double group Dj h [|| can 
be built from three rotations through angle tt about the 
coordinate axes k—x, y, z, called the signature operators, 



Rk =e 



—iirlk 



(1) 



to which one adds the inversion operator P and the time- 
reversal 



(2) 



where Ik = ^2n=i Jh^ * s ^ ne total angular momentum op- 
erator, j'[ and 0t are the angular momenta and the 
Pauli spin matrices for the particle number n, respec- 
tively, and K is the complex-conjugation operator in the 
coordinate representation. 

Following the convention introduced in p], with ro- 
man symbols, like U =R k or T, we denote operators act- 
ing in the Fock space H=H © Hi © ... © Ha ® .... 
Moreover, in order to help the reader in distinguish- 
ing between properties of these operators when they act 
in even, H+=H Q © H 2 © ... © Ha=2 P © ■ • ■, or odd, 
H-=Hi © H3 © ... © Ha=2 P +i © • • ., fermion spaces, we 
denote the former ones with bold symbols, and the lat- 
ter ones with script symbols, i.e., we formally split the 
Fock-space operators U=tj+U into two parts according 
to their domains. 

It follows |T^,[l6|,|(| that T>2 h is an Abelian group of 
16 elements, which contains: the identity E, inversion 
P, time-reversal T, their product P T =PT, three signa- 
tures Rfc, three simplexes S/c=PRfc, three T-signatures 
Rf =f R fc , and three T-simplexes S£=T S fc , i.e., 



where all these operators act in even-fermion-numbcr 
space TL+. 

Similarly, the Fock-space operators U, as well as the 
odd-fermion-number operators IA, form the group D™ 
which is a non- Abelian group of 32 elements. Apart from 
the 16 operators enumerated for D^, it contains their 
partners obtained by multiplying every one of them by 
the operators E or £, respectively. These operators can 
be identified with the rotation operators through angle 
2tt about an arbitrary axis. The partner operators are 
denoted by replacing the hats with bars, i.e., the group 
of operators acting in H- reads 



d td 



{£,P,f,P T ,1Z k ,S k ,1Zl,S?, 
£,r,f,V T ,TZ k ,S k ,TZl,S?}. 



(4) 



The complete Dj h and D™ multiplication tables have 
been given and discussed in Ref. ||, and will not be re- 
peated here. We only recall a few properties of the D™ 
group that are essential in the following analysis, namely, 



Kt = Si = T 2 =£, 
IR[) 2 =(SZY='P 2 =£, 



(5a) 
(5b) 



for k = x,y, z, 



Tz k ni = s k 3t = nlnj = sis? = n m , (6) 

for (k, l,m) being an even permutation of (x, y, z), and 

1i k iii = s k Si = 1i T k nJ = s?s? = H m , (7) 

for (k,l,m) being an odd permutation of (x,y,z). 

The multiplication table of Dj h is obtained by replac- 
ing £ and lZ m by E and R m , respectively, and using all 
bold symbols in Eqs. (§)-(0)- 

Obviously, a product of conserved symmetries is a con- 
served symmetry, and consequently, the conserved sym- 
metries form groups that are subgroups of D^ h or D™ . 
Therefore, in order to analyze various physically mean- 
ingful subsets of the conserved Dj h or D™ operators, 
we should first consider the subgroup structure of these 
groups. 

Suppose that in a given physical problem, the mean- 
field states obey the symmetries of a given subgroup 
rather than those of the whole Dj h or D™ groups. In 
such a case that subgroup contains the maximal set of op- 
erators representing the symmetry of the problem, i.e., all 
D^jj od D™ operators which do not belong to such sub- 
group are the broken symmetries. From the view-point of 
physics, we are more interested in the symmetries which 
are broken (which is related to interesting dynamical cor- 
relations), than in those which are conserved. It then 
follows that the physically interesting information will 
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be attached to the operators that do not belong to the 
subgroup studied, but do belong to D^ h or D™; those 
latter ones do not necessarily form a group. 

First we consider the single group Dj h , because (i) it 
is a smaller and simpler group than D™, and (ii) the 
operator £ which makes the difference between the single 
and the double group is always a conserved symmetry. 

The analysis below is based on identifying sets of the 
so-called subgroup generators, i.e, operators from which 
the whole given subgroup can be obtained by their suc- 
cessive multiplications. Choices of generators are, of 
course, non-unique, and hence in each case we discuss 
and enumerate all the available possibilities. 



A classification of all the 65 non-trivial subgroups of 
D^ h (we do not include trivial subgroups {E} and D T 



2h 



itself) is presented in Table Q. Every subgroup is assigned 
to a certain type, and described by a symbol given in the 
first column of the Table. The types are defined accord- 
ing to: (i) the number of generators in the subgroup (1, 
2, or 3), (ii) the number of Cartesian axes involved in the 
subgroup (0, I, or III standing for 0, 1, or 3), and (iii) the 
number of signature operators in the subgroup (A, B, or 
D standing for 0, 1, or 3). 

The classification is based on two important charac- 
teristics of each subgroup. As shown in Ref. ||, every 
conserved symmetry, labeled by one of the Cartesian di- 
rections x, y, or 2, induces a specific symmetry of local 
densities, related to this particular direction. Therefore, 
the number of Cartesian axes involved in the subgroup 
gives us the number of symmetries of local densities in- 
duced by the given subgroup. In addition, the number 
of signature operators illustrates the way in which the 
given subgroup is located with respect to the standard 
D2 subgroup, which is composed of the three signatures. 



classified according to the types described in the text. The 
second column gives the generators. The third column gives 
numbers of different subgroups irrespective of names of Carte- 
sian axes, and the fourth column gives the total number of 
subgroups in each type. 



Type 



Generators 



A. Subgroups of D 2 h 

Since the square of every element of D J h is proportional 
to the identity operator E, we have fifteen two-element, 
one-generator subgroups, each of them composed of the 
identity and one of the other Dj h operators. We denote 
these subgroups by {Gi}, where Gi is the generic symbol 
corresponding to one of the non- identity elements of D^- 
Obviously, only one choice of the generator is possible for 
every of the two-element subgroups. 

Similarly, group Dj h has 35 different four-element sub- 
groups, which can be called the two-generator subgroups, 
and are denoted by symbols {Gi,G 2 }- pertaining to 
their generators. The two-generator subgroups contain, 
in addition to Gi and G2, also the identity E and the 
product G1G2. Since this product is also one of the 
Dzii operators, we have in each of the four-element, two- 
generator subgroups three possibilities to select the gen- 
erators. 

Finally, there are 15 different eight-element, three- 
generator subgroups of D^, denoted by {Gi,G 2 ,G 3 }. 
Each of these subgroups contains the identity E, the three 
generators, three products of pairs of generators, and the 
product of all three generators. Hence, to choose the 
generators we may first pick any pair out of seven non- 
identity elements (21 possibilities), and next pick any 
other subgroup element, except the product of the first 
two, (4 possibilities). Since the order in which we pick 
the generators is irrelevant, one has altogether 28 pos- 
sibilities of choosing the three generators in each of the 
eight-element, three-generator subgroups of Dj h . 

In the same way one can calculate that there is 168 
different choices of the four generators of the whole Dj h 
group; one of them is, e.g., the set {T, P, IL X , R y }. This 
illustrates the degree of arbitrariness in implementing cal- 
culations for which the whole group Dj h is conserved. 
Similar freedom, although to a lesser degree, is available sch emes in even-fermion systems, 
when conserving any of the subgroups of Dj h . Of course, 
the freedom of choosing generators cannot influence the 
final results, however, it allows using different quantum 
numbers, phase conventions, and structure of matrix el- 
ements, as discussed in Sec. III. 



Generic Total 



1-0a 


{T}, {P}, {P J } 


3 


3 


Ha 


{ill}, {s k }, {sD 


3 


9 


1-Ib 


{B*} 


1 


3 


Total number of one-generator subgroups: 


7 


15 


2-0a 


{T,P} 


1 


1 


2-Ia 


{S fe ,f}, {Sj,P}, {Rl,P T } 


3 


9 


2-I s 


{R fc ,f}, {R fc ,P}, {R fe ,P T } 
{Si, S m } 


3 


9 


2-IIIa 


1 


6 


2-IIIb 


{R;,R m }, {R;,S m }, {R;,S m } 


3 


9 


2-IIId 


{R;, R m } 


1 


1 


Total number of two-generator subgroups: 


12 


35 


3-Is 


{Rfc.T.P} 


1 


3 


3-IIIb 


{R;,S m ,T}, {R;,S^,P}, {R;,R^,P T } 


3 


9 


3-IIIu 


{R;,R m ,T}, {Ri,R m ,P}, {R,R m ,P T }3 


3 


Total number of three-generator subgroups: 


7 


15 




Total number of subgroups: 


26 


65 



Classification of the subgroups of the single group 
allows us to discuss conserved and broken-symmetry 



3 



B. Subgroups of D 2 h 

In order to discuss the conserved and broken symme- 
tries in odd-fermion systems, we now proceed to the dis- 
cussion of the subgroups of the double group D™. In 
fact, the classification of Table | can now be repeated 
almost without change. Indeed, whenever a given Dj h 
subgroup contains the time reversal T, signature R^, or 
simplex at least one of those, the corresponding sub- 
group of D™ contains T, TZk, or Sk, and it automati- 
cally becomes a doubled D™ subgroup, with exactly the 
same generators. This is so, because in the D™ group 
the squares of the time reversal, signature, and simplex 
operators are equal to £ , Eq. (|5a|), and hence whenever 
one of these operators is present in the subgroup, it gener- 
ates the appropriate double subgroup of the double group 
D™ ■ On the other hand, when none of these generators 
are present in a given T)J h subgroup, this subgroup be- 
comes the subgroup of D™ without doubling. 

Therefore, all the Dj h subgroups listed in Table | are 
simultaneously subgroups of D™, provided the genera- 
tors denoted with bold symbols are replaced by the corre- 
sponding script generators. Most of the D™ subgroups 
have twice more elements than the corresponding Dj h 
subgroups, with few exceptions: subgroups {V}, {Ti^}, 
{S^}, V} do not double, and contain the same num- 
ber of elements as the corresponding subgroups of DA- 



TABLE II. Same as in Table g but for additional sub- 
groups of the double D™ group. 



Type Generators Generic 


Total 


2-0 a: {T,£} 1 
2-U: {&%,£}, {£%,£} 2 


1 

6 


Total number of two-generator subgroups: 3 


7 


3-Ia: {St,P,£} 1 


3 


Total number of three-generator subgroups: 1 


3 



Total number of subgroups: 4 10 



resentations M, and hence, from the view-point of the 
symmetry breaking it is irrelevant whether or not a given 
subgroup has spinor representations. 

Apart from the phase relations of electromagnetic mo- 
ments H (note that the standard multipole operators are 
defined by singling out the z axis), the three Cartesian 
directions are, of course, entirely equivalent. Therefore, 
even though changing names of axes leads to different 
subgroups of Dj h or D™, they are identical from the 
point of view of physically important features. In Ta- 
bles | and 0, we give in the third columns the numbers 
of generic subgroups, i.e., those which are different irre- 
spective of names of axes, and in the fourth column - the 
total numbers of subgroups of each type. Index k always 
denotes one of the axes, i.e., k can be equal to x, y, or 
z, while indices I and m, l^m, denote one of the three 
pairs of different axes. 

Subgroups in types "0" do not depend on the Carte- 
sian axes and, therefore, for them the numbers of generic 
subgroups are equal to the total numbers of subgroups. 
Those in types "I" have one generic form each, and three 
forms in total, depending on which Cartesian axis is cho- 
sen. Finally, for subgroups in types "III" , the total num- 
bers of subgroups can be the same, three times larger, 
or six times larger (2-IIIa) than the numbers of generic 
subgroups. 

In practical applications, conservation of different Dj h 
or D™ subgroups may require considering either only 
the generic subgroup, or all the subgroups with changed 
names of axes. For example, if one considers a triaxially 
deformed system with 0°<7<60°, the lengths of prin- 
cipal axes a y <a x <a z define the orientation of the nu- 
cleus. Then, conserved Dj h or D™ subgroups with dif- 
ferent names of axes may lead to different physical con- 
sequences. On the other hand, it can be advantageous to 
consider only one generic subgroup, with a fixed orien- 
tation, and allow for various orientations of the physical 
system by extending allowed values of 7 deformation be- 
yond the standard first sector 0°<7<60°. 



In addition, these few exceptional subgroups can be 
doubled explicitly by adding Z to the set of generators. 
For completeness, these additional subgroups of D™ are 
enumerated in Table 0. However, the physical contents 
of the additional, and of the corresponding not-doubled 
subgroups from Table ||, are the same. For example, 
they lead to exactly the same symmetry properties of 
the density matrices ||. The difference between them 
consists in the fact that the latter ones have no irreps 
in the spinor space, whereas the former ones have one- 
dimensional irreps with spinor bases (cf. p5[). However, 
from the view-point of the symmetry breaking, they lead 
to exactly the same schemes, and thus the additional sub- 
groups shown in Table || can be called trivial. Note, that 
the single-particle operators (e.g., the mean-field Hamil- 
tonian) are classified according to one-dimensional rep- 



III. SINGLE-PARTICLE BASES AND MATRIX 
STRUCTURE OF THE SINGLE-PARTICLE 
HERMITIAN OPERATORS 

Throughout this section we restrict our analysis to 
hermitian single-particle operators, and we study their 
matrix elements in the single-particle space. Therefore, 
we are here concerned with the odd number of particles 
(one), and hence we have to consider the double group 
D™. As discussed in [||, the D™ operators are either 
linear or antilinear, and they can have squares equal to 
either unity or minus unity, as summarized in Table II of 
Ref. ||. This gives us four categories of operators with 
markedly different properties, which here will be used for 
different purposes. 

In the discussion which follows, we assume that the 
single-particle basis is composed of pairs of time-reversed 
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states, and in addition we assume that the spatial wave 
functions are real, i.e., not affected by the time-reversal, 



T|nC) = Cl™ -0. 



(8) 



where <^=±1 represents the intrinsic-spin degree of free- 
dom, and n represents the set of quantum numbers cor- 
responding to space coordinates. In particular, for the 
harmonic-oscillator (HO) basis, n=(n x , n y , n z ) are the 
numbers of quanta in three Cartesian directions. As- 
sumption (j^) does not preclude whether or not the time- 
reversal is a conserved operator; it only defines the prop- 
erty of the single-particle basis in which the dynamic 
problem is to be solved. (In principle, the discussion be- 
low can be mutatis mutandis repeated with X replaced 
by V T , however, the use of the time-reversal operator is 
more appropriate in practical applications.) 

From now on we also assume that the basis is ordered 
in such a way that its first half corresponds to the (=+1 
states, and the second half is composed of their time- 
reversed (=— 1 partners. In fact, we are entirely free to 
choose states in the first half of the basis (£=+1), and 
then Eq. (0) defines those which belong to the second 
half (C=— 1). In such basis, the single-particle matrix el- 
ements corresponding to an arbitrary hermitian operator 
O have the form 



O = 



A Y 
yt B 



(9) 



where A and B are hermitian matrices, Y is arbitrary, 
and all submatrices are, in general, complex. 



A. Single-particle bases for conserved D 2h operators 

We may now separately consider several cases corre- 
sponding to different subgroups of conserved D™ op- 
erators, Table |, and to the four different categories of 
operators (linear or ant ilinear, and her mitian or antiher- 

we consider cases 



IIIA1 



[II A 4 



mitian). In subsections 
of various D™ ge nerators being separately conserved, 
and in subsections III A5 - HI A 8 , cases of pairs of gen- 
erators being simultaneously conserved. Three -generator 
subgroups are briefly discussed in Sec. Ill A 9. 



1 . Time-reversal 

Let us first consider operators O which are either even 
(invariant) or odd (antiinvariant) with respect to the 
time-reversal: 

-Pdf=e T 6, e T = ±l. (10) 

From Eqs. (||) and (|l0|) one gets 

(nC\G\n'C)=e T CC'(n ~C\0\n' -(')*■ (11) 



It follows, that the matrix corresponding to O reads 

° = ( -etY* e?A* ) > < 12 ) 

where A is hermitian, and Y is antisymmetric or sym- 
metric for £t=+1 and 6t=— 1, respectively. No block- 
diagonal structure appears, nevertheless, only two in- 
stead of four submatrices, A and Y, complex in general, 
have to be calculated. 



2. T -signature or T '-simplex 

As is well known [[l7],[|, f° r eacn of the six antilin- 
ear D™ operators, i.e., for T-signature TZT or T-simplex 
S£, k—x,y 7 z, which have squares equal to unity, (Z 2 =£, 
where Z denotes one of them), one can construct a basis 
composed of eigenvectors of Z with eigenvalue equal to 
1, 



Z\nQ = K). 



(13) 



Moreover, since every operator Z commutes with the 
time-reversal X, such basis can always be chosen so as 
to fulfill condition rt§) at the same time. Table [II lists 



examples of such bases, constructed for the HO states 
s 2 =±i). A similar construction is possible for 
any other single-particle basis, and has the explicit form 
shown in Table III provided the space and spin degrees of 



freedom are separated. Note that any linear combination 
of states |nf=+l) and \nC=— 1), with real coefficients, 
is another valid eigenstate of Z with eigenvalue 1. 
For operators even or odd with respect to Z, 



Z^OZ = e z 6, 



e Z = ±l, 



one then has 



(nC|<V O =e z (n(\6\n'Cy 



(14) 



(15) 



Hence, in bases fulfilling Eqs. (||) and (|l^), matrix O is 
purely real (ez=+l) or purely imaginary [ez=— 1). This 
gives matrix O in the form (tilde stands for the transpo- 
sition) 



A Y 
Y B 



or 



A' Y' 
-Y' —B' 



(16) 



for 6z=+l and ez——l, respectively, where all subma- 
trices are real, A and B are symmetric, A' and B' are 
antisymmetric, and Y and Y' are arbitrary. Note that 
in order to diagonalize O one only needs to diagonalize 
a real matrix with unrestricted eigenvalues (for e z =+l), 
or an imaginary matrix with pairs of opposite non-zero 
eigenvalues (for cz=— 1). 
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TABLE III. Examples of eigenstates \n x n y n z f)J of the 
T-signature, , or T-simplex, <Sj operators, k=x, y, or 
2, with eig envalue 1 [cf. Eqs. (§) and @], determined for 
the harmonic oscillator states I s 2 =±A). Symbols 

(N x , N y , N z ) refer to (n x ,n y ,n z ) for Sj operators, and to 
(Tiy+n z , n x +n z , n x +n y ) for TZ^ operators. 



k C \n x n y n z ()l 

x +1 (+i) JVa; exp(-if ) In^nj,^, s z =+|) 

a; -1 (—i) N * ex P(+*f ) \n x n y n z , s z =-\) 

y +1 (+i) JV « +1 In^n*, 

2 +1 4g (\n x n y n z , s z — \) +i(-l) Nz \n x n y n z , s z = 

2 —1 L — \\n x n y n z ,s z — ^) — ~ \n x n y n z , s z — 



real eigenvalues with no additional restrictions, or (ii) for 
€x=— 1, one hermitian matrix Y'Y, again with complex 
submatrices in twice smaller dimension, which gives pairs 
of real non-zero eigenvalues with opposite signs. 

TABLE IV. Eigenstates | C)fe of the signature or 

simplex operators, TZt or Sk for k=x, y, or z, [cf. Eqs. 
([]) and dig)], determined for the harmonic oscillator states 
\n x n y n z ,s z —±^). Symbols (N x , N y , N z ) refer to (n x ,n y ,n z ) 
for 5fc operators, and to (n y +n z , n x +n z , n x +n y ) for 7^.*, oper- 
ators. Phases of eigenstates are fixed so as to fulfill condition 



•|>) 



k c 



x 
x 



y +1 

y -i 

2 +1 
2 -1 



-±= (\n x n y n z ,s z 

(-i) N °> n 

[\nxn y nz,s z 



\n x n y n z C) k 

=5) - (-1)™ X \n x n y n z 



W) 



k)+ {-1) N,C \n x n y n z ,s z =-±}) 



i^i (Inan^ne, s»=i) l)^" |n^n y n z , s a =— ~) ) 
■ (\n x n y n z ,s z = ±) +i(-l) N * \n x n y n z , s z =-|) ) 



+r* exp(-if ) In.n^n., s z =+i(-l)^ +1 > 



exp(+if ) \n x n v n z ,s z - 



-§(-ir +1 > 



3. Signature or simplex 



Let us now consider operator which is even or odd 
with respect to one of the six linear D™ operators, sig- 
natures 1Z k or simplexes Sk, k=x,y,z, which have squares 
equal to minus unity, (X 2 =£=— E, where X denotes one 
of them), i.e., 



X^6x = e x 6, e x = ±l. 



(17) 



Since every operator X commutes with the time-reversal 
T, one can always choose a basis in which Eq. (|sj) holds, 
and 



X\n()=t(\nO, 



(18) 



where again £=±1, so r=£i=±i is the signature (for 
X=lZk) or simplex s—Qi—±i (for X—Su) quantum num- 
ber. Table |y| lists such bases constructed for the HO 
states \n x n y n z ,s z =±7j). A similar construction is pos- 
sible for any other single-particle basis. Note that one 
can arbitrarily change the phases of states \n £=+1), and 
still f ulfill Eq s. ( 3) and (18); we shall use this freedom in 
Sees. [Ill A 8| andttllBl. 



From Eqs. (U7f) and (18) one gets: 



(n(\6\n'C) =e x (C(nt\d\n'C), 
so the matrix O has the form 



A 
B 



Y 
Ft o 



(19) 



(20) 



for ex—+i and ex=— lj respectively, with A and B her- 
mitian, and Y arbitrary complex matrix. Note that in 
order to diagonalize O one only needs to diagonalize: (i) 
for ex=+l, two hermitian matrices with complex, in gen- 
eral, submatrices in twice smaller dimension, which gives 



Com paring result (|20| ) with that obtained in Sec. 
[II A 2 , one sees that the antilinear symmetries allow for 



using real matrices, while linear symmetries give special 
block-diagonal forms for complex matrices. 



4- Parity 

Standard simplification always occurs for operators 
which are even with respect to the parity operator V, 



v^ov = 6. 



(21) 



All matrices and submatrices introduced above or be- 
low acquire a block-diagonal form, provided the single- 
particle bases consist of states with well defined parity 
(such as, e.g., bases listed in Tables [II or IV). There- 



fore, apart from Sec. Ill A 7 , we do not separately discuss 
cases when the parity is one of the conserved D™ oper- 
ators, and we note that the effect of the parity conserva- 
tion can be easily included on top of any other symmetry 
conditions. 



5. Time-reversal, and T-signature or T -simplex 

For operators 6, for which both Eq. @ and (|lj) hold, 
matrix O in Eq. (12) can be additionally simplified, and 
reads 



A Y 
-ctY ctA 



A' 



Y' 



e T Y> -e T A' 



(22) 



G 



for ez=+l and ez——l, respectively, where all subma- 
trices are real, A is symmetric, A' is antisymmetric, 
Y=-e T Y, and Y'=-e T Y'. 

In particular, with t T —ez—+'\., the matrix from Eq. 
(E2) reduces to 



O 



A Y 

-Y A 



(23) 



where A is symmetric, Y antisymmetric, and both are 
real. In order to diagonalize such matrix, one can con- 
sider a smaller problem, by constructing a complex ma- 
trix Oc=A—iY which has the size twice smaller than the 
original matrix O. After diagonalizing Oq, and separat- 
ing real and imaginary parts of its complex eigenvectors, 



Oc{u + iv) = u(u + iv), 



one gets two degenerate real eigenvectors of O: 



(24) 



and 



If O is the mean-field Hamiltonian, such a 



form of eigenvectors simplifies expressions for densities. 



6. Time-reversal, and signature or simplex 

For operators O, for which both Eq. ( |l0| ) and (|l7]) hold, 
matrix O in Eq. (^) can be additionally simplified, and 
reads 





) 


or 


( ° 




o 


e T A* 


I -e T Y* 


I) 



(25) 



for ex=+l and ex=— 1, respectively, with A hermitian, 
and Y antisymmetric (ey=+l) or symmetric (£t=~ !)• 
Of course, this case is identical to that described in 



Sec. Ill A 5, because whenever the time-reversal, and sig- 
nature or simplex are conserved, the corresponding T- 
signature or T-simplex are also conserved, and ex—ez^T- 
Therefore, we may now use two different bases, and ob- 
tain two different forms of the matrix O, ( p2] ) or (25), 
which lead either to real, or to block-diagonal matrices. 
Note that in order to diagonalize matrix O for ex=+l, 
one has to diagonalize only its hermitian submatrix A, 
which has dimension twice smaller than O, similarly as 
in Eq. A 



7. Parity, and signature or simplex 

In the D™ group, the possibility of having at ones dis- 
posal two different quantum numbers simultaneously is 
very limited. Indeed in D™ one nas only three pairs of 
commuting linear operators, namely, (TZk,V) for k=x, y, 
or z. For each such pair, the corresponding simplex op- 
erator Sk is also conserved, but it does not give any addi- 
tional quantum number. Only one generic two-generator 



subgroup, {71k, V}, see 2-I# in Table |, allows, therefore, 
for two quantum numbers. Similarly, only three generic 
three-generator subgroups allow for two quantum num- 
bers, namely, (i) {TZk, 7~, V}, which allows only for sta- 
tionary solutions, (ii) {lZi,7Z m ,V}, which does not allow 
for non-zero average values of the angular-momentum, 
and (iii) {TZi,S^ l ,'P}, which is the only two-quantum- 
numbcr subgroup which allows for rotating mean-field 
states. Needless to say, this latter case is most often 
used in cranking calculations to date, see Sec. Ill C . 



8. T -signature or T-simplex, and signature or simplex 

Let us now consider operator O which is even or odd 
with resp ect to o ne of the six antilinear, Z 2 —£, operators 
(see Sec. [II A 2), and simultaneously even or odd with 
respect t o one o f the six linear, X 2 =£=— £, operators 
(see Sec. [II A 3). In such a case, simplification of the 
single-particle basis is possible only for pairs of Z and 
X operators which correspond to two different Cartesian 
directions. Indeed, focusing our attention on signatures, 
7Z T commutes with TZk , and therefore (being antilinear) 
it flips the Tlk signature quantum number. Therefore, 
the eigenstates of TZk cannot be eigenstates of TZ T . We 
may then only w ork either in the basis of eigenstates of 
7? T 



Sec. Ill A 3 



Sec. [II A 2, or in the basis of eigenstates of TZu 



On the other hand, for Ij^k TZf anticom- 
mutes with TZk, and therefore, it conserves the TZk signa- 
ture quantum number. Hence, the eigenstates of TZk can 
be rendered the eigenstates of TZf by a suitable choice of 
phases. 

It is easy to check that after multiplying eigenstates 
listed in Table [V by the following phase factors: 



$ ifc =exp{iCf(AT i + l)}, fork<l, (26a) 

S lfc =exp{tC!+*C§(JVl+iV* + l)}, forZ<fc, (26b) 

one obtains the basis states, 

\nOik = $ik\n() k , (27) 

which simultaneously fulfill Eqs. (§h, (13), and ([l8|), i.e., 



f\nC) lk =C\n -Oik, 
Zi\nC,)ik = \n()ik, 
X k \nC,)ik = «C|nC)ife. 



(28a) 
(28b) 
(28c) 



Here, Zi stands for TZj or Sf, and Xk stands for 
TZk or 4- In Eqs. ©, symbols (N x ,N y ,N z ) re- 
fer to (n x ,n y ,n z ) for the or 5^ operators, and to 
(n y +n z , n x +n z , n x +n y ) for the TZj or TZk operators. 
Moreover, a circular ordering of Cartesian directions is 
assumed, i.e., x<y<z<x, in order to define conditions 
l<k and k<l. 
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For operators 6 even or odd simultaneously with re- 
spect to Zi and Xk, see Eqs. ( |l4|) and (^7|), bases denned 
by Eq. ( p7| ) allow for a very simple forms of matrices O. 
Combining conditions ( ^6|) and ( ^p| ) one obtains block 
diagonal and real matrix elements, e.g. for ez=+l, 



A 
B 



or 



F 

y o 



(29) 



for ex=+l and ex=~ 1, respectively, with A and i? real 
symmetric, and Y arbitrary real matrix. 



9. Three-generator su* 



Apart from the unique case of the whole D™ group 
being conserved, which amounts to conserving its four 
generators, we also have 15 different three-generator 
subgroups (Table |l|), which when conserved, may lead 
to physically different mean-field solutions. Conserved 
three-generator subgroups are exceptional in that they 
do not lead to further simplifications of the matrix ele- 
ments of mean-field Hamiltonians. 



This is so, because cases enumerated in Sec. Ill A 1 
III A8| exhaust different possibilities of using conserved 
DzIF operators to simplify the structure of operators by 
suitable choices of the single-particle bases. Indeed, the 
type "III" subgroups of D™ , Table |, which involve oper- 
ators for three different Cartesian axes, do not induce any 
new simplifications. The signature or sim plex op erators 
for different axes (the X operators of Sec. [ill A 3| ) do not 
commute, and hence cannot give independent quantum 
numbers of single-particle states. Similarly, T-signature 
or T-si mplex o perators for different axes (the Z operators 
of Sec. Ill A 2| ) do not commute either, and hence cannot 
simultaneously define phases of single-particle states. 

One should stress, however, that even if a given con- 
served symmetry does not allow for any further simplifi- 
cation of the matrix elements of a mean-field Hamiltonian 
(like each third generator of a three-generator subgroup) , 
its conservation or its non-conservation may induce en- 
tirely different solutions of the mean-field problem. 



B. Phase conventions 



In Sec. Ill A we have shown how one can simplify the 



matrix elements of operators by using a given phase con- 
vention, i.e, by fixing phases of single-particle basis states 
in a given way. Whenever an antilinear D™ operator is 
conserved, one can always construct a phase convention 
for which the matrix elements of the mean-field Hamil- 
tonian are real numbers. However, from technical point 
of view, it can be more advantageous to fix the phase 
convention in yet another way. Indeed, whenever the cal- 
culation of matrix elements is more time consuming than 
the diagonalization of the Hamiltonian matrix, one may 
use the phase convention to facilitate the former task, 



at the expense of diagonalizing complex matrices. More- 
over, such a strategy allows for keeping the simplicity 
of performing the former task even in cases when there 
is no antilinear conserved symmetry available, and when 
one has to diagonalize complex matrices anyhow. In the 
present section we show constructions of phase conven- 
tions which facilitate calculations of the space-spin ma- 
trix elements. 

Representation (||) , which separates space and spin de- 
grees of freedom, is convenient in applications pertaining 
to deformed single-particle states, as those discussed in 
the present study. This is because, each hermitian oper- 
ator can be represented as a sum of four components of 
the form 



where acts in the coordinate space, and <7 M arc the 
Pauli matrices acting in the spin space, with <to defined 
as the 2x2 identity matrix. Then, the matrix elements 
of can be factorized into space and spin parts 

(nfl^VO = (n\6M\n') ■ (CMC), (31) 

and the spin part can be computed once for all. Usu- 
ally many of the spin matrix elements (CI^mIC) vanish, 
thus making it unnecessary to calculate the correspond- 
ing coordinate-space matrix elements {n\0 ( f'\n') . 

Matrix elements of operators can be made purely 
real or purely imaginary if phases of single-particle basis 
states are chosen in such a way that, for one Cartesian 
direction l=x, y, or z, one has 



A* = 0,1, 2, 3, 



(30) 



/Q|nC) = |nC), 



(32) 



where IC?=£ is the antilinear spin operator defined by 



Ki = fi&i = ia{I ' . 



(33) 



Indeed, for time-even (ey=+l) or time-odd (ct=— 1) op- 
erators one obtains that 

t\6^fc x = e^iP&^t = e^erO^, e M , - ±1, (34) 

where coefficients e M ; are given in Table |v|. Using 
Eqs. (p2[) and (p4|), one gets for matrix elements of 



(nCl^VO = Wr(n CI^VO* (35a) 
(n(\6M\n'C) = e^CC(n -C|0 (m V ~C), (35b) 



where ( 35a ) tells us whic h elements are real, and which 
are imaginary, while (35b) gives the matrix elements, e.g., 
for 1 expressed through those for (=+!• 

As is usual for antilinear operators, there is a lot of 
freedom in finding bases (32) of eigenstates of /Q. We 
can use this freedom to fulfill other useful conditions. 
For example, since AC/ and IZk anticommute for l^k, one 
can find bases ( |32"| ) which are at the same time the eigen- 
states of signature or simplex operators. In fact, phases 
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of eigenstates listed in Table [V , has been chosen in such 
a way that, 



operator (p3|), and calculations were performed within 



lCi\n() k = \n() k , for k < I, 



(36) 



where again the circular ordering of Cartesian directions, 
x<y<z<x, is assumed to define k<l. 

TABLE V. Antilinear spin operators fti, which can be 
used to fix convenient phase conventions leading to conditions 
(§)• 



/ 


tl 


eoi 


eii 


til 


(31 


X 




+ 1 


+ 1 


-1 


-1 


y 


Tio y =K, 


+1 


-1 


+1 


-1 


z 


Tia z =—ia x IC 


+ 1 


-1 


-1 


+ 1 



C. Examples of previous cranking approaches 



As argued in Sec. [II A 7 , there are good reasons to 
use in cranking calculations the subgroup {TZi 1 Sj n ,'P} 
(Table |) of conserved D™ symmetries. This generic 
three-generator subgroup appears in three space orienta- 
tions, i.e., for l=x, y, or z, and each of these possibilities 
was employed in one of the HF(B) or phenomenological- 
mean-field cranking analyses to date. 

In particular, traditionally the x axis was chosen as 
the direction of the cranking angular momentum, see 
e.g. Ref. JllJl , and therefore, the standard Goodman ba- 
sis Q corresponds to the l—x subgroup, with phases of 
single-particle states (and quasiparticle states, for that 
matter) fixed by using the operator. Then, by 
dropping the parity operator from the symmetry group 
{1Z X , V}, most octupole-cranking calculations were 
performed within the 2-IILa subgroup {S X ,S^} of Table 

B 

Another choice was made in the HO-basis [|J and 
coordinate-space [ft2]|ll| HF(B) calculations, where the 
z axis was used as the cranking axis. Such choice was 
motivated by the standard representation of spinors, 
that are eigenstates of a z , and hence the l=z subgroup 
{TZ Z , Sy , V} was employed. In these approaches, phases 

of single-particle states were fixed by using the Sy opera- 
tor, and the parity-broken calculations were done within 
the {S z ,Sy} subgroup. 

Finallyin the recent Cartesian HO-basis HF approach 
of Ref. fl4| , the code HFODD was constructed for the 
conserved l=y subgroup {7Z y , Sj, V}, and the y direction 
was used for the cranking axis. The choice of this sym- 
metry, and the resulting choice of the y cranking axis, 
was motivated by the fact that it allows for using real 
electric multipole moments, cf. Ref. ||. Phases of single- 
particle states were in Ref. || fixed by using the JC Z 



the basis of the S y eigenstates, Table gVJ. The HFODD 
code allows for calculations with one symmetry plane, 
and this is done within the {S y } conserved symmetry 
group of Table |. The code can also optionally perform 
the two-symmetry-plane cranking calculations for the 2- 
III A subgroups {S y ,S^} and {S y ,Sj}. 



IV. CONCLUSIONS 

We have analyzed the "far end" of the symmetry 
breaking chain, namely, symmetries of mean-field nu- 
clear states, which range from time-even, parity-even, 
signature conserving states, (nevertheless breaking the 
rotational and axial symmetry), to those which do not 
conserve any symmetries at all. We have shown that in- 
termediate cases, between such two extremes, correspond 
to conserved subgroups of the Dj h or D™ point symme- 
try groups. A classification of all the subgroups has been 
proposed, and we have shown that there are 26 different 
non-trivial symmetry-breaking schemes, when names of 
Cartesian axes are irrelevant, and 65 different non-trivial 
symmetry-breaking schemes when names of axes are dis- 
tinguished in the intrinsic frame of reference. 

Consequences of conserving individual D™ symme- 
tries have been enumerated for the construction of single- 
particle bases in which mean-field operators may have 
special simplified forms. We point out that the same 
forms of the mean-field Hamiltonian may correspond to 
different conserved symmetries, and hence to different 
physical consequences for observables obtained in the 
mean-field methods. We have also analyzed and com- 
pared various options for defining phase conventions of 
single-particle basis states. 
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